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We study anomalies in the Coulomb blockade spectrum of a quantum dot formed in a silicon 
nanowire. These anomalies are attributed to electrostatic interaction with charge traps in the 
device. A simple model reproduces these anomalies accurately and we show how the capacitance 
matrices of the traps can be obtained from the shape of the anomalies. From these capacitance 
matrices we deduce that the traps are located near or inside the wire. Based on the occurrence of 
the anomalies in wires with different doping levels we infer that most of the traps are arsenic dopant 
states. In some cases the anomalies are accompanied by a random telegraph signal which allows 
time resolved monitoring of the occupation of the trap. The spin of the trap states is determined 
via the Zeeman shift. 
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I. INTRODUCTION 

Single electron charges or spins are very appealing as 
logic bits, either as ultimate classical bits or quantum bits 
if coherence is usedA*£. To read such bits, either quan- 
tum point contacts or single electron transistors (SETs) 
are used. SETs have indeed been used as very sensitive 
electrometers for the (time-averaged) charge on a second 
quantum dot for over a decade nowAi More recently 
the radio-frequency SET technique 5 was used to moni- 
tor the charge on the second dot or to measure a current 
by electron count ing&£i2i2. This allows to measure lower 
currents than standard measurements and gives access to 
the full counting statistics of the current. 

Such experiments are difficult because any device that 
involves detection of few or single electron charges is 
subject to the dynamics of surrounding charge traps. 11 
This is particularly critical for metallic SETs™ For SETs 
based on the very mature silicon CMOS technology the 
control of this offset charges seems to be better. 13 

These charge traps are quantum dots whose presence 
or properties are not controlled. Typically they consist 
of defects on atomic scale. Their sizes are therefore much 
smaller than what is possible for lithographic quantum 
dots. If their positions, although being random, can be 
limited to some zone, the charge traps are not necessarily 
a nuisance but can useful. An example are flash memo- 
ries where the trend is to replace the lithographic floating 
gate by grown silicon nanocrystals inside the gate oxide. 
They are grown in a layer and have all the same distance 
from channel and gate electrode. Their dynamics are 
therefore very similar Another example are dopants 
in semiconductors. Efforts are made to control their in- 
dividual position in a silicon crystal. 16 Indeed, besides 
the location, their properties are very uniform and solid 
state quantum bits based on dopants in a silicon crys- 



tal -individually addressed by gates and contacts- were 
proposed as solid state quantum bit a 17 i 18 i 19 . Silicon is 
interesting as host material because the spin relaxation 
time can be very long 20 compared to GaAs. The detec- 
tion of spins of individual traps in a silicon field effect 
transistor has been recently reported using random tele- 
graph noiso 21 i 22 . However, in this experiment the traps 
seem to be in the oxide rather than in the silicon. 

In this work, we use nanowire-based silicon transistors 
operated as SETs at low temperature to detect the lo- 
cation, spin and occupation number of individual charge 
traps, which we attribute to As dopant states. They are 
capacitively coupled to the SET. Therefore they induce 
anomalies in the otherwise very regular periodic oscilla- 
tions of the drain-source conductance G versus gate volt- 
age V g . We compare the data with simulations obtained 
after solving the master equation for the network formed 
by the main dot and the charge trap. 

Not only the static time-averaged current is analyzed 
but also the switching noise which appears near the de- 
generacy point in gate voltage where the trap occupation 
number fluctuates. 

Finally, a magnetic field was applied in order to probe 
the spin polarization of the traps via their Zeeman shifts. 
As expected from simple considerations 23 , we observed a 
majority of singly occupied traps. 



II. SAMPLES AND SETUP 

Samples are produced on 200 mm silicon on insula- 
tor (SOI) wafers with 400 nm buried oxide and a boron 
substrate doping of 10 15 cm -3 . The SOI film is locally 
thinned down to approximately 20 nm and a 30 nm wide 
and 200 nm long nanowire is etched from it. A 40 nm 
long polysilicon control gate is deposited in the middle 
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FIG. 1: (color online) Sample layout and electrical model. 
The insert shows a top view of the sample before back-end 
process obtained in a scanning electron microscope. The main 
image shows a transmission electron micrograph (TEM) of a 
type B sample along the silicon nanowire (black, the wire 
shown here is thinner than in the samples used for measure- 
ments). Light gray regions are silicon oxide. The darker re- 
gion in the center is the polysilicon gate with SisN4 spacers on 
both sides of it. Below, a schematic energy diagram is drawn. 
The reduced doping level below the spacers and the gate elec- 
trode creates a potential barrier, in the middle of which a well 
is created by a positive gate voltage. Conductance through 
the barriers separating the well from source and drain occurs 
by tunneling through a chain of well connected dopants (plot- 
ted in the right barrier)^ In more isolated dopants (plotted in 
the left barrier) the number of charges is well quantified. Such 
traps are the main concern of this paper. Their interaction 
with the quantum well is mainly electrostatic. We describe it 
with the lumped network superimposed to the TEM. 



of the wire (see Fig. QJ. There are two layouts. Type A: 
The wires are uniformly doped with As, above 10 19 cm -3 . 
The gate oxide is 4nm. Type B: The wires are first uni- 
formly doped at a lower level (As, 10 18 cm -3 ), then, after 
deposition of the gate electrode and 50 nm- wide S13N4 
spacers on both sides of it, a second implantation pro- 
cess increases the doping to approximately 4 • 10 19 cm -3 
in the uncovered regions while the doping level stays low 
near the gate. In this layout the gate oxide is 10 or 24 nm 
thick, with a 2 or 4 nm thermal oxide and 8 or 20 nm de- 
posited oxide. Most measurements are made on type B 
samples, and we use type A samples mainly for compar- 
ison. 

The measurements were performed in a dilution refrig- 
erator with an electronic base temperature of approxi- 
mately 150 mK. We used a standard 2- wire low frequency 
lock-in technique with low enough voltage excitation to 
stay in the linear regime and a room temperature cur- 
rent amplifier (gain 100 MSI). For time resolved measure- 
ments a DC bias voltage was applied and current mea- 
sured with a 10 current amplifier (bandwidth 10 kHz) 
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FIG. 2: (color online) Drain-source conductance versus gate 
voltage for 3 different samples. All samples have the same 
width and gate length but the sample in the upper panel is of 
type A with a 4 nm gate oxide, while the samples in the lower 
panels are of type B. The one of the middle panel has 10 nm 
gate oxide, the one of the lower panel 24 nm. The smooth 
field-effect characteristics at room temperature (black lines) 
are replaced by Coulomb blockade oscillations at base temper- 
ature (blue curves). The period is determined by the surface 
area of the nanowire/gate overlap. The Coulomb blockade 
oscillations in the upper panel are irregular compared to the 
ones in the lower panel where only some rare anomalies per- 
turb the very regular spectrum. These anomalous regions 
with reduced contrast and fluctuating peak spacing are high- 
lighted with circles. The anomalies marked with bold circles 
are studied in detail in this work. 



followed by a 33 kHz AD conversion. For spin sensitive 
measurements a superconducting magnet was used to ap- 
ply an in-plane magnetic field up to 16 T. 



III. DATA 

Figure [2 shows typical G(V g ) plots. At room temper- 
ature our samples behave as classical (albeit not opti- 
mized) n-channel MOSFETs. Below approximately 20 K 
they turn into single electron transistors with regularly 
spaced Coulomb blockade resonances. The period Vc — 
-gg of these oscillations (e is the absolute value of the 
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electron charge) is determined by the gate capacitance 
C g , which in turn can be estimated from the gate/wire 
overlap and the gate oxide thickness. For the sample 
with 4nm (A), 10 nm (B), 24 nm (B) gate oxide the peak 
spacing is respectively 14 mV ± 4 mV, 10.3 mV ± 0.5 mV, 
15.3 mV ± 0.8 mV. This corresponds to a gate capaci- 
tance of 11 aF, 15.5 aF, 10.5 aF. For the type A samples 
the gate capacitance is in good agreement with the simple 
planar capacitance estimation. For the type B samples 
where the gate oxide thickness is of the same order as the 
dimensions of the wire, the 3-dimensional geometry has 
to be taken into account. The gate capacitance of the 
type B samples is increased with respect to the type A 
samples because the flanks of the wire play a more im- 
portant role. A 3-dimensional numerical solution obtains 
a good agreement with the measured capacitances. 

The peak spacing statistics has already been measured 
and compared to theory^ Here we focus on anomalous 
regions where the conductance contrast is markedly re- 
duced and a phase shift of the Coulomb blockade oscilla- 
tions occurs. This results in tails in the Gaussian peak- 
spacing distribution. Such perturbations to the periodic 
pattern are marked with circles in Fig. O In the type B 
samples with low doping, these perturbations occur only 
rarely (we observe typically 3 to 5 per sample). In the 
unperturbed regions, the height of the Coulomb blockade 
peaks shows long-range correlations. In the type A sam- 
ples with high doping level the perturbations are more 
frequent and let the whole spectrum look irregular, (see 
top panel of Fig. |2J) This suggests that the perturbations 
are related to the doping. 

In the measured stability diagram, i.e. the 2D plot of 
conductance versus gate and bias voltages, the pertur- 
bations are even more visible (see Fig. |3J). In the per- 
turbed regions additional teeth appear in the Coulomb 
diamonds. 

We develop a simple model based on a trap state lo- 
cated in the vicinity of the quantum dot, and compare 
the simulation with the experimental data. 



IV. MODEL 

The quantum dot formed by the gate electrode in the 
middle of the wire is separated from the source and drain 
reservoirs by a piece of silicon wire containing only a 
few tens (type B) or hundreds (type A) of dopants. In 
the type A samples these access regions extend from the 
border of the gate electrode to the regions where the wire 
widens (see Fig. [Q and its resistance becomes negligible. 
In the type B samples only the zones below the spacers 
contribute significantly to the access resistance and the 
highly doped parts of the wire can be considered as part 
of the reservoirs. 

Electrons pass through these access regions by trans- 
port via the dopant states. 42 . As the dopants are dis- 
tributed randomly and the coupling between them de- 
pends exponentially on their distance, this coupling is 
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FIG. 3: 2D-plots of the measured drain-source conductance 
versus gate and drain voltages in an unperturbed, very peri- 
odic gate voltage range (upper panel), and in an anomalous re- 
gion where a charge trap is observed (lower panel). White ar- 
eas correspond to Coulomb blockaded regions (no detectable 
current). The lines inside the conducting regions are not the 
excited states of the dot (the spacing being too high and al- 
most identical for all resonances) . We attribute these lines to 
additional conduction channels on the drain side that open at 
higher bias (chains of well connected dopants lying somewhat 
higher in energy than the drain Fermi level). Compared to 
the measurements in Fig. |5| the anomalous region has shifted 
by 50 mV in gate voltage after thermal cycling between base 
and room temperature. We do not observe such shifts as long 
as the sample is kept cold. 

distributed over a wide range. Transport therefore takes 
place mainly through a percolation path formed by well 
connected dopants 2 ^ while other dopant states are only 
weakly connected and their occupation is a good quan- 
tum number (see Fig. [IJ. We attribute the anomalies 
in the Coulomb blockade spectrum to the electrostatic 
interaction of the quantum dot with such a charge trap 
formed by an isolated dopant site. 

We model this with the lumped network shown in 
Fig , n] Similar models have been considered in Refs. 
l26l and l2?l A small trap (t) is capacitively coupled to 
source (s), gate (g) and to the main dot (m). We note 
d = Cf + Cf + Cf and X % = CfV s + C?V d + CfV g 
(i = m, t) After some calculation, the electrostatic energy 
of the two dot system can be expressed as a function of 
the charges Q m on the main dot and and Q t in the trap. 

M(Q m ,Q t ) T(Q t ) 

where C c is the capacitive coupling between dot and trap, 

A = c^fc> c = Cm + and x = Xm + 
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For a small trap (C m < C t ) these renormalizations are 
weak: C « C m and X w X m . The problem is symmetric 
under exchange of main dot and trap even though the 
expression in Eq. n °t. W is plotted in the top panel 
of Fig. m 

We focus on the structure of the Coulomb blockade 
conductance fixed by Eq. JJ) and not on the exact value 
on the conductance plateaus. Therefore we choose as 
simple as possible the following parameters which are 
necessary for the simulation but do not affect the struc- 
ture of the conductance diagram. 

We suppose all transmission coefficients to be constant, 
the ones connecting the main dot to source and drain be- 
ing 1000 times higher than the ones connecting the trap 
to the main dot and source or drain. Electrons can there- 
fore be added or removed from the trap, but their con- 
tribution to the total current through the device is negli- 
gible. This contrasts with models of stochastic Coulomb 
blockade 28 or in-series quantum dots22^& where the cur- 
rent has to pass through both dots. 

In terms of kinetic energy, we describe the main dot 
as metallic (negligible single-particle level spacing A, i.e. 
A <C kT) and we consider only one non-degenerate en- 
ergy level for the trap. In source and drain we suppose 
a uniform density of states. We assume fast relaxation 
of kinetic energy inside the dot and the reservoirs, i.e. 
thermal distributions in the electrodes and the main dot, 
even for nonzero bias voltage. With these assumptions, 
the transition rates of an electron in the main dot to 
the source or drain reservoirs or from the reservoirs to 
the dot are proportional to the auto-convolution of the 
Fermi function. The transition rates from or towards the 
trap are directly proportional to the Fermi function™ 

The statistical probability for each state (Q m5 Qt) of 
the system can now be calculated by solving numerically 
the master equation and gives access to the mean current 
through the system. 

Results of such a numerical study are presented in 
Fig. 21 The middle panel shows the mean occupation 
of the trap. On a large scale, the trap becomes occupied 
with increasing gate voltage. In the central region of the 
figure however, whenever an electron is added onto the 
main dot, the electron in the trap is repelled and only 
later it is re-attracted by the gate electrode. Inversely, 
the trap charge repels the charges on the main dot and 
the Coulomb blockade structure of the main dot is shifted 
to higher gate voltage when the trap is occupied (see 
lower panel). The two Coulomb blockade structures for 
unoccupied and occupied trap are respectively indicated 
by dotted and dashed lines in the middle and lower panel 
of Fig. H 

This explanation is illustrated in terms of energy in 
the top panel of Fig. 21 which shows the energies for 
the different charge configurations. The crossings of the 
blue (black) parabolas give the positions of the Coulomb 
blockade peaks for empty (occupied) trap. The shift be- 
tween the crossings of the black parabolas with respect 
to the crossings of the blue parabolas and the shift of 
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FIG. 4: (color online) Numerical study of a trap coupled 
to the source and to the main quantum dot, as sketched in 
FigQ] Parameters: effective temperature: T = 1 K; main dot: 
Cg = 60 % , d = C s m = 70 f ; trap: Cf = 0.045 f , C s t = 2 f , 
Cf = 0, Co = 1 v • The trap can either be empty or charged 
with one electron. The upper panel shows the energy for the 
different charge states at zero bias in function of gate volt- 
age. The blue parabolas are for empty trap, the black ones 
for occupied trap. The thick blue and black lines indicate the 
ground state of the main dot for respectively empty and oc- 
cupied trap. The middle panel shows the self-consistent mean 
occupation number of the trap, the lower panel the resulting 
conductance through the dot. The effect of the charge trap 
is to shift the Coulomb blockade diamonds of the main dot 
depending on the charge in the trap. The dotted (dashed) 
lines indicate the position of the diamonds when the trap is 
empty (occupied). This result is in very good agreement with 
the experimental data shown in the lower panel of Fig. |3] 

the dashed lines with respect to the dotted lines are 
due to the term /3 t Qt in M(Q m ,Q t ). Knowing that one 
Coulomb blockade oscillation corresponds to a change of 
e in /3 t Qt + X, the shift due to AQ t = e is 

SV g = frVc (2) 

where Vq is the Coulomb blockade peak spacing of the 
main dot. 

We will now determine the width of the anomaly in 
the Coulomb blockade spectrum at low bias voltage. It 
is given by the gate voltage range where the occupation of 
the trap oscillates at zero bias. In the top panel of Fig. 0] 
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this is the zone between the first and the last crossing 
of the thick black line and the thick blue line. First we 
calculate AM, the difference of the ground state energies 
for empty and occupied trap arising from the term M 
m Eq.ffl Then we calculate the change in gate voltage 
necessary for T(—e) — T(0) to exceed this difference. 

AM reaches its extreme values when for one state of 
the trap main dot is at a degeneracy point (the kinks in 

the thick lines), where M = ^c - * ^ or ^ ne °^ ner state 
of the trap the main dot is a fraction A of a Coulomb 
blockade period away from the degeneracy point and 

M = e2(1/ 2 2 ^ t)2 - The extrema of AM are therefore 

±^A(i-A)- 

The gate voltage dependence of term T is given by 
at = ^3%(T(-e)-T(0)) = jgfa. Note that a t is 
the long-range gate voltage lever arm of the trap over 
several Coulomb blockade oscillations, where the charge 
of the main dot has to be considered as relaxed with the 
source and drain Fermi levels. T(—e) — T(0) has to pass 

form +^/?t(l - A) to - fcA(l - A) m order to toggle 
the trap definitively. The width AV g of the anomaly is 

therefore given by — ea t AV g = — 2^A(1 — A) or 

n = A(l-A)— (3) 

n is the number of anomalous periods and a m = ^ with 
C g = + AC? the gate voltage lever arm of the main 
dot. 

We have identified at = c -tq and A = c+c as 
parameters that determine the structure of the trap sig- 
nature. Both do not depend on the absolute value of 
the trap's capacitances. Indeed, if one allows only or 

1 electron in the trap, the absolute value of the trap 
capacitances enters the problem only indirectly by mod- 
ifying slightly the capacitance matrix of the main dot 
and cannot be obtained in the limit of a small trap. Our 
model contains therefore only 2 effective parameters for 
the trap instead of 3 (Cf , C t s , C c ). All 3 parameters of 
the trap are only significant if the trap can accommodate 

2 or more electrons. In this case the spacing between the 
anomalies gives access to the absolute values of the trap's 
capacitances 

Figure illustrates the relation between the trap's ca- 
pacitance matrix and its signature. If teeth of constant 
width for all anomalous resonances are visible at the pos- 
itive slope of the Coulomb blockade diamonds, the trap 
is on the source side of the dot. If they are visible at 
the negative slope, the trap is on the drain side. The 
width of the teeth depends on A? the width of the anoma- 
lous region essentially on a t (for A close to - where the 
anomalies are well visible). 

V. POSITION AND NATURE OF THE TRAPS 

As an illustration, from the lower panel of Fig. [3| we 
infer a t « 0.015 and A ~ 0-3. These are the actual 
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FIG. 5: Calculated trap signatures for different sets of pa- 
rameters, (a) The trap is close to the source, (b) The cou- 
pling to the gate electrode is reduced by a factor of 5. The 
signature becomes wider, (c) The coupling to the source is 
reduced, the coupling to the dot increased, (d) The trap is 
placed on the drain side of the dot instead of the source side, 
(e) The coupling to the dot is reduced, the coupling to the 
drain increased. 



parameters that have been chosen for the simulation in 
Fig. 21 and the lower panels of Fig. [3| and Fig. 2] are in- 
deed very similar. As for all impurities we observed, a t 
is small. This is what we expect for a trap inside the 
silicon wire. The coupling to the gate electrode is much 
weaker than the coupling to the main dot or the source 
electrode because the dielectric constant of the oxide bar- 
rier (esio 2 = 4) is much smaller than that of bare sili- 
con (esi = 12), which in addition is enhanced near the 
insulator-metal transition. 32 

Positions of the trap outside the wire can be ruled out. 
Traps located deep inside the oxide can be excluded be- 
cause their transmissions would be too weak to observe 
statistical mixing of occupied and unoccupied trap states 
during our acquisition time below 1 s. Similar devices in- 
cluding intentional silicon nanocrystals at the interface 
between thermal oxide and deposited oxide have been 
studied in views of memory applications™^ The mea- 
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sured lifetime of charges in the nanocrystals exceeds 1 s 
by orders of magnitude already at room temperature and 
at low temperature gate voltages of about 5 V have to be 
applied in order to toggle the charge in the nanocrys- 
tals. The traps must therefore be inside the Si wire or at 
its interface with the oxide. But the interface traps are 
unlikely. We estimate their density to be smaller than 
10 11 cm -2 , corresponding to a few units per sample. As 
they are distributed throughout the entire band gap it 
is very unlikely to observe several of them in the small 
energy window a t (V g ma * - V g min ) « 30 meV that we scan 
in our measurement. The most likely traps are therefore 
defects in the silicon wire or As donor states. Given the 
volume of the access regions under the spacers and the 
doping level iVo, there are around 70 donor states under 
the spacers in devices of type B. We estimate the width 



of the impurity band to be 
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150 meV. 23 One 



should therefore expect around 15 dopants in the energy 
window. Typically we record 3 to 5 anomalies. Indeed 
we do not expect to observe anomalies for all dopants 
because the charge on well connected dopant sites is not 
quantified and, according to our model, dopants very 
close to the dot (f3 t « 1) or to the reservoir ((3 t « 0) 
produce very small anomalies. 



In the type A samples the doping level in the access 
regions is more than 10 times higher than in the type B 
samples. The whole Coulomb blockade spectrum should 
therefore be anomalous. Indeed, the spectrum is much 
less regular (see Fig. than for the type B samples, es- 
pecially for low gate voltage, but we cannot distinguish 
signatures as clear as in the type B samples. This is 
consistent because in the type A samples the mean dis- 
tance between impurities is less than 3 nm and they are 
too well connected for the charge on them to be well 
quantified. In other words, the wire is very close to the 
insulator-metal transition. Our doping level is in fact 
already higher than the bulk critical As concentration 
N c = 8.6 • 10 18 cm~ 3 M^ 



We have deduced that the observed traps lie inside the 
wire. The position of the trap along the wire can also be 
determined. Traps on the source side and the drain side 
of the dot can be distinguished (see Fig. [SJ the teeth of 
constant width SV g appear on the positive slope of the 
diamonds in case of a trap on the source side of the dot 
and on the negative slope in case of a trap on the drain 
side) and the parameter (3 t gives the ratio between the 
capacitances towards the main dot and the source (or 
drain) electrode. As the dielectric constant of the wire 
is much higher than the surrounding silicon oxide, this 
ratio can be translated linearly to a position in direction 
of the wire. In the example of Fig. with flt « 0.3 we 
would expect the impurity to be located | on the way 
from the dot (border of the gate electrode) to the source 
reservoir (source side border of the spacer). 
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FIG. 6: (color online) Analysis of a trap signature with 
switching. Sample of type B with 24 nm gate oxide. The 
width of the wire is 80 nm instead of 30 nm. (a) A RTS 
trace taken at ]/ g = 500 mV, Va = — 6mV. Light gray 
trace: raw data. Black trace: data after compensation of 
the time constant of the current amplifier. Blue line: fitted 
signal, (b) Histograms of the times passed in the weak cur- 
rent state (occupied trap, black) and the —1 nA state (empty 
trap, blue). The time constants (averages of these times) are 
0.31ms and 0.62 ms. The corresponding exponential distri- 
butions (straight lines) fit well the histograms, (c) Current 
histogram at Vd = — 6mV. The nonzero density between the 
two current levels is due to the finite rise time. The current 
for unoccupied trap is always higher than for occupied trap, 
(d) Time constants of the empty and occupied levels in func- 
tion of gate voltage and occupation number of the trap. 



VI. TIME-RESOLVED OCCUPATION NUMBER 

In the preceding sections we assumed charge traps with 
changing mean occupation number to explain anomalies 
in the mean conductance through a Coulomb blockaded 
quantum dot. Yet the measurements of the mean current 
have not allowed us to measure the occupation number of 
the trap directly. But the currents through the main dot 
differ for empty and occupied trap because the position 
of the Coulomb blockade resonances is shifted, and at 
the anomalies where the occupation number of the trap 
is different from and 1, the fluctuations of the occu- 
pation number should create a random telegraph signal 
(RTS) 21 i 22 i 36 i 37 in the current through the main dot. 

Indeed, we frequently observe strongly increased cur- 
rent noise near the anomalies, especially at low gate volt- 
age. However, at most anomalies at higher gate voltage 
we do not observe a clear increase of the noise level. This 
indicates that changes in the trap state occur at frequen- 
cies far beyond 10 kHz, the bandwidth of our measure- 
ment. Indeed, for charge traps formed by dopants we 
would expect the transmission rates of the trap to be of 
the same order as for the main dot, where the transmis- 
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sion occurs also through dopant states. The excess noise 
of the trap should therefore be comparable to the shot 
noise of the quantum dot, which does not emerge from 
the noise floor of the current amplifier. But much smaller 
transmissions should also be possible as the dopants are 
distributed randomly and the transmission rate depends 
exponentially on the distance between them. In some 
cases we observe clear RTS with time constants larger 
than 100 us. An example is given in Fig.UJa). The distri- 
bution of the times spent in the two states follows the ex- 
ponential distribution expected for a RTS (see Fig.JSfb)). 

The color plot of the current distribution in Fig. EJc) 
shows the evolution of the the two current levels (dark 
lines with high probability) with gate voltage. Above 
380 mV the two levels are very different. This difference 
is most likely due to electrostatic interaction of the trap 
and the current path through the barrier: depending on 
the state of the trap, the dopants through which the main 
part of the current flows are well or poorly aligned in en- 
ergy. The fact that the current levels never cross simpli- 
fies greatly the assignment of the high and low current 
levels to the states of the trap. The high current trace 
being most likely at low gate voltage and the low cur- 
rent trace being most likely at high gate voltage allows 
to attribute the high current to empty trap and the low 
current to occupied trap. 

The time constants of the empty and occupied state 
are plotted in Fig. Efd). Consistently with panel (c), 
the time constant for the empty trap decreases with gate 
voltage while the time constant for the occupied trap 
increases. Superimposed with this slow change there are 
oscillations with a period of 12 mV, the peak spacing of 
the main dot. This oscillation is even more prominent 
in the mean occupation number given by r ° ccu f ed . 

^"occupied T^empty 

As explained in section HVl for the case of low bias, this 
oscillation is due to the discrete charge on the main dot 
which cycles the trap several times between empty and 
occupied state. It is not observed in RTS in larger devices 
without Coulomb blockade^. At high bias (Vd > only 
an oscillation of the occupation probability remains of 
this cycling. This can be seen in Fig. Cfb) and (d) where 
the occupation probability for different bias voltages is 
compared with simulation. As in Fig. the oscillations 
in Fig. [7f b) and (d) are aligned parallel to the negative 
slopes of the Coulomb blockade diamonds indicating that 
the trap is on the source side of the dot. 

RTS (i.e. current through the trap) only occurs when 
the trap is in the bias window. For large gate and bias 
voltage excursions where the charging energy of the main 
dot is negligible, the main dot can be considered as part 
of the drain reservoir. The zone where the trap is in 
the bias window is then delimited by slopes -^h and — 
(indicated by straight lines in Fig.0), just as for a single 
quantum dot. These slopes give a more straightforward 
access to the parameters a t and f3 t . 

The mean occupation of the trap is higher for positive 
drain voltage than for negative drain voltage indicating a 
higher transmission rate of the trap towards source than 




FIG. 7: (color online) Comparison of measured occupation 
number and simulation. Same trap as in Fig. El (a) Mean 
differential conductance obtained by numerical derivation of 
the mean current, (b) Occupation of the trap obtained from 
the duty cycle of the RTS signal. Regions where no clear 
RTS could be detected are left white, (c) and (d) Simulation 
with the following parameters: main dot: = 80 |r, Cm = 
60 f , C m = 100 f ; trap: C? = 0.08 f , C t s = 0.6 f , C t d = 0, 
C Q — \^j. In units of the drain-dot barrier transmission, the 
source-dot barrier transmission is 10 for empty trap and 
for occupied trap, the source-trap barrier transmission 
and the trap-dot barrier transmission o^q- 



towards the main dot. 

In Fig.[7fc) and (d) we try to reproduce panels (a) and 
(b). For this simulation we reduce by a factor of 100 the 
transmission of the source barrier of the main dot when 
the trap is occupied. This reproduces the lines of reduced 
differential conductance at positive drain voltage (com- 
pare Fig.[7Ja) and (c)). In the simulation the oscillations 
of the trap occupation decay more rapidly with bias volt- 
age than in the measurement. This could be related to 
our approximation of a thermal distribution of kinetic 
energies in the main dot, which is certainly not accurate 
at high bias voltage. 

Charge traps are generally believed to be not only re- 
sponsible for RTS noise but also for 1/f noise in SETsii 
and decoherence^. These interpretations imply a large 
number of traps with small influence on the device (in 
our model f3 t ~ 0). Such traps could be dopants in the 
reservoirs or the substrate. 
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VII. SPIN 

The spin of the trap state leads via the Zeeman energy 
under magnetic field to a gate-voltage shift of the trap 
signature of: 



| I I 111 | I I II I | ■! I I I M I I | I 



eCXt ~dB = 9flB 



(4) 



fiB is the Bohr magneton, AS Z the change in spin quan- 
tum number of the trap state in direction of the mag- 
netic field when an electron is added to the trap. It can 
take the values ±~. If there are already electrons in 
the trap higher changes are also possible, but they imply 
spin flips and such processes are therefore expected to be 
very slow^ The Lande factor g for impurities in Si and 
SiC>2 has been measured by electron spin resonance4& 
The observed renormalizations are beyond the precision 
of our measurements, therefore we take g = 2. The gate- 
voltage lever-arm of the trap states at is very weak as we 
have shown above. The Zeeman shifts should therefore 
be strong. 

Indeed, the magnetic field clearly shifts the trap signa- 
ture in Fig. [SI to lower gate voltage. In order to identify 
the shift as the Zeeman effect, we compare it quantita- 
tively with the prediction of our model. The shift of the 
resonances due to the trap is half the peak spacing, so 
/?t = \ (see Eq. 0). The lever arm for the main dot is for 
this gate voltage a m = 0.26 and the width of the trap sig- 
nature varies from 2.5 periods without magnetic field to 
1.5 periods at 16 T. This implies a gate- voltage lever arm 
for the trap of a t = 0.026 ••• 0.043 (see Eq. Q) which 
we interpolate as a linear function of magnetic field. The 
dotted line in Fig. [HI is obtained if we put this lever-arm 
and S z = — \ in Eq. It is in very good agreement with 
the measured shift and confirms our model. The increase 
of the lever arm with magnetic field could be explained as 
follows. In the access regions the nanowire is close to the 
metal-insulator transition and the dopant states strongly 
increase the dielectric constant 34 . Under magnetic field 
they shrink 23 , reducing the localization length and the 
dielectric constant in the wire. Therefore the coupling 
towards the main dot and the reservoir decreases while 
the gate capacitance dominated by the oxide capacitance 
remains unaffected. 

We observe such Zeeman shifts in the majority of our 
samples. In most cases the trap signature shifts to lower 
gate voltage as in Fig. [HI This is what we expect for 
isolated traps occupied with one electron. When a trap 
state is occupied with a second electron it has to oc- 
cupy the energetically less favorable state whose energy 
is increased by the Zeeman effect. This leads to a shift 
towards higher gate voltage under magnetic field. Al- 
though isolated As-donor sites in Si can only be occupied 
by one 43 electron due to Coulomb repulsion, clusters of 
two donors could contain two or more electrons 4i For not 
too high doping levels clusters should however be rare. 
Accordingly we observe much less shifts to higher than 
to lower gate voltage. In devices based on similar technol- 




i i i i i i i 
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FIG. 8: (color online) Shift of a trap signature with magnetic 
field. The dotted line indicates the Zeeman shift expected for 
a trap state being occupied by a first electron. It depends on 
the gate- voltage lever-arm which in turn is determined by the 
width of the signature. This prediction of the Zeeman shift 
follows exactly the observed shift. 



ogy Xiao et al. observed that all shifts occurred to higher 
gate voltag o 21 i 22 indicating doubly occupied traps. With 
precise measurements of the Lande factor they located 
the traps inside the oxide. This difference also supports 
that the traps in our device are not located in the oxide 
but inside the silicon wire. 



VIII. PERSPECTIVES 

Dopant states in silicon could provide very scalable 
solid state quantum bits, based on charge, electron spin 
or nuclear spin. But it is still very difficult to control 
their position individually. On the other hand, with 
several gate electrodes one could imagine to select suit- 
able dopants out of large number of randomly distributed 
dopants. In this context we have presented how the ca- 
pacitance matrix of charge traps near a small silicon sin- 
gle electron transistor can be determined and we showed 
how the gate- voltage dependence of the occupation is re- 
lated to the spin of the trap state and that the charge 
in these traps can be read out. These charge traps are 
attributed to arsenic dopant states. At a doping level of 
10 18 cm -3 we observe several well isolated dopant states 
per device as well as percolation paths of well connected 
dopants linking the main quantum dot to the reservoirs. 
In similar geometries with multiple gate electrodes the 
coupling between the dopants could be tuned by chang- 
ing their alignment in energy with the well connected 
dopants. Such randomly distributed dopants are proba- 
bly more suited for electron spin quantum bits than for 
charge quantum bits where two dopant sites with small 
distance are necessary. In this perspective we are work- 
ing on measurement of the coherence time of the elec- 
tron spin in the observed traps. Together with the excel- 
lent stability in time as well as its full compatibility with 
CMOS technology our system could be a good basis for 
scalable quantum bits. 
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